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ABSTRACT 

We provide a new analytic blastwavc solution which generalizes the Blandford-McKee solution to 
arbitrary ejecta masses and Lorentz factors. Until recently relativistic supernovae have been discovered 
only through their association with long duration Gamma Ray Bursts (GRB). The blastwaves of such 
explosions are well described by the Blandford-McKee (in the ultra relativistic regime) and Scdov- 
Taylor (in the non-relativistic regime) solutions during their afterglows, as the ejecta mass is negligible 
in comparison to the swept up mass. The recent discovery of the relativistic supernova SN 2009bb, 
without a detected GRB, opens up the possibility of highly baryon loaded mildly relativistic outflows 
which remains in nearly free expansion phase during the radio afterglow. In this work, we consider 
a massive, relativistic shell, launched by a Central Engine Driven EXplosion (CEDEX), decelerating 
adiabatically due to its collision with the pre-cxplosion circumstcllar wind profile of the progenitor. 
We compute the synchrotron emission from relativistic electrons in the shock amplified magnetic field. 
This models the radio emission from the circumstcllar interaction of a CEDEX. We show that this 
model explains the observed radio evolution of the prototypical SN 2009bb and demonstrate that SN 
2009bb had a highly baryon loaded, mildly relativistic outflow. We discuss the effect of baryon loading 
on the dynamics and observational manifestations of a CEDEX. In particular, our predicted angular 
size of SN 2009bb is consistent with VLBI upper limits on day 85, but is presently resolvable on VLBI 
angular scales, since the relativistic ejecta is still in the nearly free expansion phase. 
Subject headings: gamma rays: bursts — supernovae: individual (SN 2009bb) — shock waves - 
radiation mechanisms: non-thermal — techniques: interferometric 
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1. INTRODUCTION 

Ultra relativistic bulk motion of matter particles in as- 
trophysical settings is i mplied most notab ly in Gamma 
Ray Bursts (GRBs) fsee lPiranl (|1999l . 12004) for reviews) . 
Afterglows from GRBs are generated from the emission 
by relativistic shocks that result from slowing down of a 
relativistic shell by the the medium surrounding the pro- 
genitor star that exploded. Similar interaction of stellar 
material (ejecta) from an exploding star with the circum- 
stcllar matter (CSM) gives rise to non-relativistic shocks 
in core collapse supernovae. 

Fluid dynamical treatment of ultra-relativistic spher- 
ical blast waves mediated by strong shock fl has been 
given by IBlandford fc McKed (fl976l Mf% . They de- 
scribe a similarity solution of an explosion of a fixed 
amount of energy in a uniform medium. This includes 
an adiabatic blast wave and an impulsive injection of en- 
ergy on a short timescale as well as an explosion where 
the total energy increases with time, suggesting that the 
blast wave has a continuous central power supply. An- 
other important model considered by them is that of 
a blast wave propagating into a spherically symmetric 
wind. On the other hand, the initial nearly free ex- 
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2 Here the strong shocks are those where the random kinetic 
energy per particle behind the shock is much greater than the un- 
shocked medium, i.e., ^2/^2 2> Pl/ni, the subscripts 2 and 1 de- 
note the shocked gas and unshocked gas respectively and p, n arc 
the pressure and number densities. 



pansion of a non-relativistic supernova blast wave, inter- 
acting wi th the sur r oundi n g circumst e llar m edium, was 
found bv IChevalierl (fl982l ): iNadezhinl (fl985l ). Once the 
blast wave sweeps up more CSM material than its own 
rest mass, the self-similar solutions of non-relativistic 
bla st waves are described in the Newtonian regime by 
the lSedovl (fl946h Ivon Neumaiml (fl96l iTavlorl (fl950[ ) so- 
lution. 

In this paper we provide an analytic solution of the 
standard m o del of relativistic hydro dynamics (see e.g. 
iPiranl 119991 iChiang fc Dermerl 1 1999T) for an adiabatic 
blastwave. Here, the exploding shell decelerates due to 
inelastic collision with an external medium. That is, we 
provide the solution for an arbitrary Lorentz factor of the 
expanding supernova shell. The need for such a solution 
which can handle a trans-rclativistic outflow is motivated 
by the discovery of SN 2009bb, a type Ibc supernova 
without a detected GRB which shows clear evidence of 
a mildly relativistic ou tflow powered by a central engine 
(ISoderberg et al.l l2niTt) . SN 2009bb-like objects (Cen- 
tral Engine Driven Explosions, hereafter CEDEX) dif- 
fer in another significant way from classical GRBs: they 
are highly baryon loaded explosions with non-negligible 
ejecta masses. Our ne w analytic blastwave solu tion 
therefore generalizes the IBlandford fc McKeel (|1976| ) re- 
sult, in particular their impulsive, adiabatic blast wave 



in a wind-like p cx 



CSM. 



The new class of relativistic supernovae without a de- 
tected GRB, e.g. SN 2009bb, relaxes two important 
and well-known constraints of the GRBs, namely the 
compactness problem and baryon contamination. Thus 
highly baryon loaded mildly relativistic outflows can re- 
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main in the nearly free expansion phase during the ra- 
dio afterglow. Therefore we consider the evolution of a 
massive relativistic shell launched by a CEDEX, experi- 
encing collisional slowdown due to interaction with the 
pre-explosion circumstellar wind of the progenitor. We 
calculate the time evolution of the radius, bulk Lorentz 
factor and thermal energy o f the decelerating blastwave . 
Our solution reduces to the IBlandford k McKed (| 1976( 1 
solution, in the ultra-relativistic, negligible ejecta mass 
limit. We also quantify the density of accelerated elec- 
trons, amplified magnetic fields and hence the radio syn- 
chrotron emission from the blastwave of a CEDEX. 

In Section [5] we outline the conditions that require 
ultra-high bulk Lorentz factors and very low ejecta 
masses in classical Long GRBs. We argue why these 
constraints are relaxed in the case of mildly relativistic 
outflows detected in SN 2009bb like events, that have no 
detected GRBs. It provides the motivation for the ana- 
lytic solution derived here. In Section [3] we develop the 
analytic solution of a relativistic blast wave launched by a 
CEDEX, for the collisional slowdown model described by 
iPiranl (j 1 9991 1 . We use this solution to show that the SN 
2009bb blast wave is substantially baryon loaded and re- 
mains in the nearly free expansion phase throughout the 
~ 1 year of observations (Figure [1]). In section @] we dis- 
cuss the relativistic blast wave energetics. We quantify 
the amount of shock accelerated electrons and magnetic 
field amplification in a CEDEX blast wave. In Section 
[5] we use this information to model the radio spectrum 
and light curve of a CEDEX in the nearly free expansion 
phase. In section [5] we provide expressions to deduce 
the blast wave parameters from the observed radio spec- 
trum. We also compare our blastwave solution to other 
known solutions relevant for relativistic blastwaves in the 
literature (Section [7|). In section [8] we discuss the impli- 
cations of baryon loading in determining the evolution 
and observational signatures of a CEDEX. In the Ap- 
pendix we provide the analytic expressions for the tem- 
poral evolution of the blast wave parameters (Scction [A"|) . 
We demon strate that our solut i on red uces to the ultra- 
relativistic lBlandford k McKed ()1976f l solution for a con- 
stant velocity wind, in the low mass, ultra-relativistic 
limit (Section [B]). In Section [C] we provide "stir-fry ex- 
pressions"0 for the radio spectrum of a CEDEX with 
known initial blast wave parameters. In Section [D] we 
invert the problem and provide handy expressions for es- 
timating the initial parameters for a CEDEX from radio 
observations. In Section |E] we provide an expression to 
compute angular size of a CEDEX from multi-band radio 
spectrum. Using the radio spectrum of SN 2009bb, we 
show that our predicted a ngular size is consisten t with 
the reported upper limits (jBietenholz et aH 12010( 1 from 
VLBI, but should be resolvable presently. 

2. RELATIVISTIC OUTFLOW IN GRBS AND CEDEX 

Relativistic supcrnovae have, until recently, been dis- 
covered only through their temporal and spatial associ- 

3 We note that a similar term "TV Dinner Equati ons" appeared 
in the Astrophysical literature in a paper on GRBs (Rhoads 1999) 
to denote equations where "numerical values for physical constants 
have been inserted, so they are ready to use without further prepa- 
ration"; in the present work, "stir-fry" expressions would allow 
users to pick and choose usable formulae derived or approximated 
in this work and quickly (re)assemble their own combination of 
useful input parameters. 



ation with long duration GRBs. An acceptable model 
for GRBs must find a way to circumvent the compact- 
ness and baryon contamination problems. The observed 
rapid temporal variability in GRBs imply a very com- 
pact source . The comp actness p r oblem pointed out by 
iRudermanl (|1975| 1 and iSchmidtl (|1978f l. indicates that 
7-ray photons of sufficiently high energies will produce 
electron-positron pairs and will not be able to come out 
due to the resulting high opacity. However, the ob- 
servation of high energy photons from GRB s was rec- 
oncile d with such theo r etical constraints by iGoodmanl 
(|1986t l and iPaczvnskil (|1986l ). allowing the high en- 
ergy photons to come out from a relativistic explosion. 
This requires initial bul k Lorentz fa ct ors of the radi- 
ating shell, 70 > 10 2 l|Piranl 1 1999]) . iShemi k Piranl 
(|1990f l pointed out that in the presence of even small 
amounts of baryons or baryon contamination, essen- 
tially the entire energy of the explosion gets locked up 
in the kinetic energy of the baryons, leaving little en- 
ergy for t he electromagnetic disp lay. This problem was 
solved by iRees k Meszarosl (119921) and independently by 
several authors dNaravan et al. r fl992l IRees k Meszarosl 
119941 : iPaczvnski k Xulll994l ). by considering the recon- 
version of the kinetic energy of the fireball into radia- 
tion, due to interaction with an external medium or via 
internal shocks. However the allowed initial mass is still 
very small (M rj 1O -6 M0), as too much baryonic mass 
will slow d own the exp losion and it will no longer be 
relativistic (jPiranl Il999) . Hence, observation of a short 
bright pulse of 7-ray photons in GRBs require a very 
small amount of mass to be ejected with a very high 
bulk Lorentz factor. 

The burst in the GRB itself results from the conversion 
of kinetic energy of ultra-relativistic particles or possi- 
bly the electromagnetic energy of a Poynting flux to ra- 
diation in an optically thin region. An inner engine is 
believed to accelerate the outflow to relativistic speeds, 
although the engine may remain hidden from direct ob- 
servations. The "afterglow" on the other hand results 
from the slowing down of a relativistic shell on the exter- 
nal medium surrounding the progenitor star. There can 
also be an additional contribution to the afterglow from 
the inner engine that powers the GRB, since the engine 
may continue to emit energy for longer duration with a 
lower intensity and may produce the earlier part of the 
afterglow, say the first day or two in GRB 970228 and 
GRB 970508 (|Piranlll999t Ikatz et al.lH^9^ l. 

In supernovae associated with GRBs, after a brief 
high energy electromagnetic display, the relativistic 
ejecta continues to power a long lived radio afterglow 
(|Kulkarni et aLlU99a ISoderberg et alJl2004D . Since the 
emergence of 7-ray photons, early in its evolution, con- 
strains the initial relativistic ejecta mass to be very small 
and the initial bulk Lorentz factor to be very large, the 
relativistic ejecta sweeps up more circumstellar material 
than its own rest mass by the time the radio afterglow 
is detected. The ev olution of the radia tive blastwave has 
been described by ICohen et al.l (|1998l ). During the ra- 
dio afterglow phase if radiative losses do not take away 
a significant fraction of the thermal energy, the blast- 
wave may be treated as adiabatic. Under these condi- 
tions, the evolution of the blastw ave is well described 
by the IBlandford k McKed (|1976f) solution if the blast- 
wave remains ultra-relativistic or by the Sedov- Taylor 
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solution if the blastwave has slowed down into the New- 
tonian regime. The interaction of GRBs with their cir- 
cumst ellar wind has been discussed by Chevali er fc Lil 
(|2000( ). The spectra and light curves of GRB after- 
glo ws have been com puted in the ultra-relativistic regime 
bv iSari et al.1 (119981) and in the Newtonian regime by 
iFrail et al.l (I2000D 

iSoderberg et al.1 (|2010D have recently discovered bright 
radio emission, associated with the type Ibc SN 2009bb, 
requiring a substantial rclativistic outflow powered by 
a central engine. The search for a 7-ray counterpart, 
spatially and temporally coincident with SN 2009bb, 
in data obtained from the all-sky Inter Planetary Net- 
work (IPN) of high energy satellites, did not detect 
a coincident GRB (ISoderberg et al.l 12010ft . The ab- 
sence of detected 7-rays from this event relaxes the con- 
straints of very high Lorentz factors and very low mass 
in the relativistic ejecta. In fact, the ra dio emitting out- 
flow in SN 2009bb is mildly relativistic (jSoderberg et al.l 
120101 ) and has a large baryonic mass coupled to it, as 
evidenced by the nearly free expansion for ~ 1 year 
(Figure [T]). Since the swept up mass is still smaller 
than the rest mass of the relativistic ejecta, the evolu- 
tion of the blastwave c an neither be described by the 
the rapidly decelerating iBlandford fe McKed (119761) so- 
lution n or the Sedov- Taylo r solut ion. iDermer fc Chiang! 
(|1998fk iChiang fc Dermerl (|1999l ) have considered the 
emission from a collisionally decelerating blast wave 
as it transitions from nearl y free expansion to 
the IBlandford fc McKed (j!976f ) self similar solution. 
iHuang et al.l (|1999l ) have discussed the subsequent tran- 
sition to a Sedov-like phase. 

Hence, the discovery of relativistic supernovae, such 
as SN 2009bb, without detected GRBs, relaxes the con- 
straints from the compactness problem and baryon con- 
tamination. This motivates the study of highly baryon 
loaded mildly relativistic outflows which can remain in 
the nearly free expansion phase during the radio after- 
glow. 

3. RELATIVISTIC BLASTWAVE SOLUTION 

We use t h e simple collisiona l mod el described by 
iPiranl (jl999l ); IChiang fc Dermerl (|1999ft where the rela- 
tivistic ejecta forms a shell which decelerates through 
infinitesimal inelastic collisions with the circumstcllar 
wind profile. The initial conditions are characterized 
by the rest frame mass Mq of the shell launched by a 
CEDEX and its initial Lorentz factor 70. In contrast 
to the self similar solutions, which describe the evolu- 
tion a way from the boundaries of the independent vari- 
ables (iBarcnbl att fc ZerDovichlll972|) and the need for 
proper normalization to get the correct total energy, our 
set of initial conditions directly fixes the total energy as 
E = ~y a M a c 2 . 

3.1. Equations of Motion 

The shell slows down by a sequence of infinitesi- 
mal inelastic collisions with the circumstcllar matter. 
The swept up circumstcllar matter is given by m(R). 
Conservation of energy and mom entum give us (see 
IChiang fc Dermedll^|plrar]ll999l ). 

dj dm 



and 



dE 



! (7 — l)dm, 



(2) 



respectively, where dE is the kinetic energy converted 
into thermal energy, that is the energy in random mo- 
tions as opposed to bulk flow, by the infinitesimal colli- 
sion. In the adiabatic case this energy is ret ained in the 
shell and we have the analytic relation (from |Piranlll999l ) 



m(R) 



= -( 7 o-l) 1/2 (7o + l) 1/2 

x /V-ir^y + ir 3 / 2 ^'- 



(3) 



At the age of 17 days, the radio luminosity of SN 
2009bb was L u ~ 5 x 10 28 ergs s" 1 Hz" 1 at 8.4 GHz 
(jSoderberg et al.lfeOlOD . This implies a radio luminosity 
of vL v ~ 4.2 x 10 38 ergs s _1 . Compared with the energy 
in rel ativistic electrons of 1.3 x 10 49 ergs (jSoderberg et al.l 
[20l0h . this gives a radiation timescale of ~ 10 years at 
age 17 days, which is very large compared to the age of 
the object. Hence, it is safe to compute the dynamics of 
the object in the adiabatic limit. We also note that since 
the available thermal energy initially grows with time as 
E oc R oc t (shown later) and the luminosity falls of as 
vL v cx t~ l (shown later) the adiabaticity condition only 
gets better with time. However, the radiation timescale 
would have been same as the age at t ~ 70 seconds, which 
is comparable to the duration of very Long Soft GRBs 
and XRFs. 

For a circumstellar medium set up by a steady wind, 
where we expect a profile with p cx r -2 , we have m(R) = 
AR where A is the mass swept up by a sphere per unit 
radial distance. A = M/v W i n d can be set up by a steady 
mass loss rate of M with a velocity of v W i n d from the 
pre-explosion CEDEX progenitor, possibly a Wolf Rayet 
star. Note that our defin ition of A i s equ ivalent to Airq 
in Equation (3 and 5) of ICheval ier (1982;). Integrating 
the right hand side then gives us 



AR _ 7V7g^l 
Mo ^fi^l 



7o- 



(4) 



For the physically relevant domain of 7 > 1 this equation 
has only one analytical root at 

7 Af + AR 

7 = ATTn— = . „, , . =FF^> ( 5 ) 



2Aj RM Q + A 2 R 2 



which gives the evolution of 7 as a function of R. T his 
is similar to Equation 8 of IChiang fc Dermerl (|1999t ) for 
a adiabatic spherical blastwave in a wind like (p oc r~ 2 ) 
CSM profile. The amount of kinetic energy converted 
into thermal energy when the shell reaches a particular 
R can be obtained by integrating Equation ([2]) after sub- 
stituting for 7 from Equation ((5|) and dm = AdR, to 
get 



E = c 2 (-M - AR 



(6) 



Ml + 2A<y RM + A 2 R 2 



3.2. Evolution in Observer's Time 
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The evolution of R and 7 can be compared with ob- 
servations once we have the time in the observer's frame 
that corresponds to the computed R and 7. For emission 
along the line of sight from a bla stwave with a constant 7 
the commonly used expression (Mcszaros fc ReeslfT997l) 
is d 

tabs = 7T~2~- (') 

27^ c 



However, [Sari (1993) has pointed out that for a deceler- 
ating ultra-relativistic blastwave the correct t b s is given 
by the differential equation 

dR 

dt bs = - 9 ■ (8) 

We substitute 7 from Equation ([5]) and integrate both 
sides to get the exact expression 

R(M + A lQ R) 

tobs ~ 2c 7 o(7oM + Ai?)- l9j 



Note that, this reduces to the iMeszaros fc Reesl ()1997l ) 
expression only in the case of nearly free expansion and 
deviates as the shell decelerates. In the rest of the work 
we use t to indicate the time i & s in the observer's frame. 
Inverting this equation and choosing the physically rele- 
vant growing branch, gives us the analytical time evolu- 
tion of the line of sight blastwave radius, as 



R = 



1 



2A l0 



x (-M + 2Acj t 



(10) 



-y/8AcM t^ + (M - 2Ac~f ty 



in the ultra-relativistic regime. This can now be substi- 
tuted into Equations J5] and [6]) to get the time evolution 
of the Lorentz factor 7 and the thermal energy E (see 
Appendix). This gives us a complete solution for the 
blastwave time evolution, parametrized by the values for 
70, Mo and A. 

3.3. Series Expansions 

Even though an analytical solution is at hand, it is 
instructive to look at the Taylor expansions in time for 
the relevant blastwave parameters of radius 

2 4(Ac 2 7o 3 (7o 2 -l)U 2 , ,v 

R = 2cy$t--± ^2 L^ + o(t 3 ), (11) 

Mo 

Lorentz factor 

^ + (t 2 ) , 



7 = 7o 



and thermal energy 



Mo 



(12) 



E = 2Ac 3 (7o - l) 7o 2 i 

2^c 4 (37o-2) 7 o 3 (7o 2 -l)t 2 
M 



(13) 



O t 3 



Equation (TT2")) immediately tells us that the blastwave of 
a CEDEX starts out in a nearly free expansion phase, 
but slows down significantly by the time tdec, when the 
first negative term in the Taylor expansion for 7 becomes 
equal to 70, where 

M 

* dec= 2(Ac7o(7 2 -l))- ( } 

This signals the end of the nearly free expansion phase 
and the solution enters a Blandford McKee like (see 



Appendix [B]) or Sedov like phase, depending upon the 
Lorentz factor at that time. 

4. BLASTWAVE ENERGETICS 

We shall now use the blastwave solution developed in 
the previous section to predict the radio evolution of a 
CEDEX with a relativistic blastwave slowing down due 
to circumstellar interaction. For the prototypical SN 
2009bb, the blastwave was only mildly relativistic at the 
time of the observed radio afterglow. In the absence of 
a significant relativistic beaming, the observer would re- 
ceive emission from the entire shell of apparent lateral 
extent Ri a t at a time t bs given by 
,, dRi a t 



/37C 



(15) 



which is valid even in the mildly relativistic regime. Be- 
cause of its mildly relativistic outflow, all expressions de- 
rived for SN 2009bb, use the above Equation rather than 
Equation [5] which is applicable in the ultra-relativistic 
case. In either case radio observations of SN 2009bb mea- 
sure essentially the transverse Ri a t not the line of sight 
R. Integrating term by term gives us the time evolution 
of the lateral radius as 



Rio 



7o 



If 



(16) 



(Ac^l^j^l)? 
Mo 



O t J 



The thermal energy available when the shell has moved 
out to a radius R is given exactly by Equation how- 
ever it is again convenient to look at its Taylor expansion, 



E 



Ac 2 { 10 - l)R 
(A 2 c 2 ( 7o 2 -l))i? 2 



(17) 



0{R A 



2M 

4.1. Electron Acceleration 

ISari et al.1 (fl99"l give the minimum Lorentz factor j m 
of the shock accelerated electrons as 

( V — 2\ m n 
7m = e e (?-— ^7. (18) 

At the mildly relativistic velocities seen in SN 2009bb, 
the peak synchrotron frequency of the lowest energy elec- 
trons are likely to be below the synchrotron self absorp- 
tion frequency. This explains the v 5 ^ 2 low frequency be- 
havior of the spectrum. Hence, considering a electron 
distribution with an energy spectrum NoE~ p dE, which 
we assume for simplicity to be extending from j m m e c 
to infinity, filling a fraction / of the spherical volume of 
radius R, we need an energy 



E, 



4/(7 m m e c 2 ) 1 -P7V 7ri? 3 
3(p-l) 



(19) 



If a fraction e e = E e /E of the available thermal energy 
goes into accelerating these electrons, then for the lead- 
ing order expansion of E in R the normalization of the 
electron distribution is given by 

AT _ 3^c 2 e e (7o - l)(7 m m e c 2 ) 2 

*° 27^ 7 (20) 

for p = 3, as inferre d from the optically th in radio spec- 
trum of SN 2009bb (jSoderberg et alj[20loh . 
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In the rest of this work we have made the simplifying 
assumption of a constant e e . The physics of shock accel- 
eration is unlikely to change during time of interest as 
relevant parameters such as the shock velocity does not 
change much. Moreover, we show later in this work that 
the light curve of a CEDEX is controlled by the com- 
petition between a decreasing synchrotron flux at high 
frequencies and a decreasing optical thickness (hence in- 
creasing flux) at low frequencies. The major light curve 
features of a shifting peak frequency and a nearly con- 
stant peak flux, arise from this effect rather than micro- 
physical processes that may change e e . 

4.2. Magnetic Field Amplification 

If a magnetic field of characteristic strength B fills the 
same volume, it needs a magnetic energy of 

E B = \B 2 fR 3 . (21) 

If a fraction es = Eb/E goes into the magnetic energy 
density, then the characteristic magnetic field is given by 



B 



c 6Ae B ( 7 o-l) 



R 



f 



(22) 



In the rest of this work, following the argument in the 
previous sub-section, we have made the simplifying as- 
sumption of a constant es- This explains the observed 
B cx R" 1 behavior (Figure [2]) seen in SN 2009bb. These 
observations strengthens the case for a nearly constant 
€b- This feature of an explosion within a p cx r~~ 2 wind 
profile is also seen in several non-relativistic radio su- 
pernovae. The normalization is given by the progenitor 
mass loss parameter, initial Lorentz factor, filling factor 
of the electrons and the efficiency with which the ther- 
mal energy is used in amplifying magnetic fields. Note 
that this phase would last only as long as the expansion 
is nearly free or t < tdec- Therefore the observed B — R 
relation in SN 2009bb argues for a tdec 1 year and 
hence a CEDEX with a highly baryon loaded outflow. 

Note that the highest energy to which a cosmic ray 
proton c an be accelerated is determ ined by the B R 
product (IHillasi 119841: IWaxmanl 120051) . We have argued 
elsewhere (jChakraborti et aLll2010f ) that the mildly rel- 
ativistic CEDEX are ideal for accelerating nuclei to the 
highest energies in sufficient volumetric energy injection 
rates and independent arrival directions to explain the 
post GZK cosmic rays. The above B cx dependence 
seen in Fig implies that a CEDEX like SN 2009bb will 
continue to accelerate cosmic rays up to a nearly con- 
stant very high energy for long times, i.e. for the entire 
time that this (nearly free expansion) phase lasts. 

5. ELECTRON SYNCHROTRON RADIO SPECTRUM 

It is expected that the radio emission from a CEDEX 
will be produced by synchrotron emission of shock ac- 
celerated electrons in the shock amplified magnetic field 
quantified in the previous section. O nce we have the 
evolution of Nq and B, we can follow IChevalierl (fl998h 
to express the Synchrotron Self Absorbed (SSA) radio 
spect rum from a CEDEX i n the optically thick regime 
using iRvbicki fc Light manl (|1979| ) as 



v ~ ce Ui 



5/2 



(23) 



where c\ , C5 and cq are constants given by iPacholczvU 
(|1970l) and D is the distance to the source. Similarly, 
the optically thin flux is given by 



-(p-l)/2 



(24) 



In mildly relativistic cases such as SN 2009bb, we may 
receive radiation from the entire disk projected on the 
sky, hence R is to be understood as Ri a t- Substituting 
for Nq and B from Equations (| 201 and [22]) and the leading 
order expansion for the projected lateral radius (Eq. I16[) , 
we have the optically thick flux as 



F u — fc 2 c 5 ( 7 o-l)(7o + W^ 



5/2^ 



(25) 



2 / 24Ae B 



1/4N 



,70/ + / 

Similarly, the optically thin flux can be expressed as 
24A 2 c 3 c 1 c 5 e B ee(la - l) 2 {j m m e c 2 ) 2 



F v 



(26) 



These equations together provide the flux density of a 
CEDEX as a function of time and frequency. 

5.1. SSA Peak Frequency 

The transition from the optically thick to optically thin 
regime happens at the peak frequency v p . At a fixed 
observation frequency this is reached in time t p . The 
condition for the SSA peak may be obtained by equating 
the optically thin and thick flux, to get 



v p t p 



223/14g5/l 



A 9 t 



cc e e e 



lof + f 



\A?(7o 2 -l) 



x(7 m m e c 2 



tt/(7o + If 



2/7 



(27) 



All quantities on the right hand side are constants or 
parameters of the problem which are uniquely fixed for 
a particular relativistic supernova. Hence, during the 
nearly free expansion phase v p cx t~ x and the peak moves 
to lower and lower frequencies as the plasma expands and 
becomes optically thin with time. This behavior is also 
seen in the radio spectra of SN 2009bb. 

The range in the values of v p t p is then weakly depen- 
dent on the initial bulk Lorentz factor 70 and the mass 
loss rate, parametrized by A. The above expression can 
be used to select the frequency and cadence of radio 
follow-ups of type Ibc supernovae for detecting CEDEXs. 

5.2. SSA Peak Flux 

The peak flux density can now be obtained by substi- 
tuting the expression for v p t p into that for F„, to get, 



F 



2 19/14 3 9/14 c 19/7 C5e 9/14 e 5/7 (A/(7o _ ^19/14 

:( 7m m e cY 0/ V/ 7 )/( C ^ 2 / 2 ). (28) 
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This peak flux is obtained by keeping only the leading 
orders, for almost free expansion, is nearly a constant for 
for t <C tdec- As the blastwave slows down, the decay of 
the peak flux can be obtained exactly from the analyti- 
cal expression for R(t). However for simplicity a Taylor 
expansion can be written down as 



F vp (t) 



= 1 - 



5 (Ac 7o 3 ) t 



Mn 



o [fi- 



lm 



The peak flux may also come down if synchrotron losses 
for the electrons are significant. This would require a self 
consistent modeling of the time dependent acceleration 
and synchrotron losses of the relativistic electrons in the 
shock amplified magnetic field of a CEDEX blastwave. 

The peak flux density F vp is then strongly dependent 
on the bulk Lorentz factor 70 and the parametrized mass 
loss rate, A. This suggests a large range in radio fluxes of 
CEDEXs and hence increased sensitivity of the Extended 
Very Large Array (EVLA) should help in detecting more 
of these enigmatic objects. 

5.3. Synchrotron Cooling 

IRvbicki fc Lightmanl ([1979D give the characteristic syn- 
chrotron frequency of an electron with Lorentz factor 
7 e 3> 1 in a magnetic field B as 



"(7e) = 77 



2 QeB 



(30) 



27rm e c 

iSari et all (|1998|) provide the critical electron Lorentz 
factor 7c, above which an electron will loose a signifi- 
cant portion of its energy within the age of the object, 

aS 67rm e c 

7c = ( 31 ) 

Substituting this into Equation 1301 we get the syn- 
chrotron cooling frequency as 

18irm e cq e 

"• = "M = -^BW- ^ 
Substituting the value of SN 2009bb mag n etic fi eld at 20 
days from Table 1 of iChakraborti et al.l (|2010f ). we get 
the cooling frequency as ~ 2.6 THz. Since the magnetic 
field decays as B cx t -1 , the cooling frequency will grow 
as v c cx t. Hence, synchrotron cooling is unimportant for 
electrons radiating near the SSA peak frequency v p ~ 7.6 
GHz, and they are indeed in the slow cooling regime 
during the timescale of interest. 

The scaling relations for v c and v p show that they 
should have been comparable at around ~ 1 day and 
at a frequency of ~ 140 GHz. Very early time (~ 1 
day) millimeter wave observations of CEDEXs are there- 
fore encouraged. They may reveal the presence of a syn- 
chrotron cooling break in the radio spectrum. This may 
be used to determine the magn etic field independent of 
the equipartition argument (see lChandra et aLll2004al) . 

6. BLASTWAVE MODEL INVERSION 

The input parameters of the model are specified by 
Mo, 70, A, e e and 6b- Under the assumption of equipar- 
tition of the thermal energy between electrons, protons 
and magnetic fields we have e e = €3 = 1/3. We then 
eliminate A between Equations ([2~T1 and |2"5)) to get 



7o 2 - 



3c 8 6 e B {D 2 F 1 



ir 8 c 9 5 e e f(j m m e c 2 )'' 



2/19 



Cl 



{Vpt p )c 



(33) 



which gives us the initial Lorentz factor 70 in terms of 
the early-time (t <C tdec) SSA peak frequency v p and the 
peak flux F vp . In the non-relativistic limit this gives us 



Zc%t B {D 2 F vp f 
ir 8 cle e f(j m m e c 2 y 



1/19 



Cl 



(34) 



which is insensitive to the equipartition parameter a = 
e e /e B as « cx a -1 / 19 . The inferred 70 has a weak depen- 
dence on 7 m and is to be solved self consistently with 
Equation (JT5J) giving us 7o ~ 1.16 for SN 2009bb. Not 
doing so will incur an error in estimating the initial bulk 
Lorentz factor. Equations (|3"31 or IMf may be used as a 
direct test of whether a radio supernova has relativistic 
eject a. 

The equipartition parameter a for a given CEDEX 
may be determined directly from radio observations if 
a synchrotron co oling break is s een in the broad band 
radio spectrum. iChandra et al.l (|2004a[) have used this 
to determine th e magnetic field in SN 1993J and its 
radial evolution (jChandra et al.ll2004b[) . independent of 
the equipartition argument. Subsequently A (essentially 
the scaled mass loss rate) may be obtained by fitting 
Equation ([22]) to the B-R data (e.g. Figure [2]). A direct 
method of estimating A would be to eliminate 70 between 
equations (J27] and [2t| to get A ~ 1.2 x 10 12 gram cm" 1 . 
For a typical Wolf Rayct wind velocity of 10 3 km s _1 this 
density profile may be set up by a constant mass loss rate 
of M ~ 1.9 X 10" 6 M© yr" 1 . This is consistent with the 
mass mass loss rate of the SN 2 09bb progenitor already 
determined by iSoderberg et all (|2010T ), 

All CEDEXs at a given distance, in their nearly free 
expansion phase, are characterized by one of the constant 
70 curves (Figure |4j and one of the constant M curves 
(again Figure |4|) . This F vp vs v p t p diagnostic plot for 
CEDE Xs is essen t ially the relativistic analogue of Fig- 
ure 4 of lChevalierl (|1998f ) where non-relativistic radio su- 
pernovae of different types are seen to lie on one of the 
constant velocity curves. For example, SN 2009bb, lies 
on a intermediate mass loss parameter A and a mildly 
relativistic 70. For the same initial 70 as SN 2009bb, an 
object encountering lower circumstellar density, would be 
significantly fainter and faster evolving, such as those oc- 
curring near the points of intersection of the solid green 
curve and the dashed red curve. These objects, possibly 
at the faint end of the luminosity function of CEDEXs, 
may be discovered by a dedicated high sensitivity search 
for relativistic outflows from nearby type Ibc supernovae 
with the EVLA. 

The ejecta mass can then be inferred from the decel- 
eration timescale tdec fitting the time evolution of the 
observed radii with the Taylor expansion (Equation I16jl 
of Riat obtained from our model. In Figure[T]we compare 
predictions from our model with the observed temporal 
evolution of the blastwave radius. The observed radii are 
consistent with a nearly free expansion for an apparent 
lateral velocity of 7/3 = 0.527 ± 0.022c. As the ejecta has 
not yet slowed down significantly, it is not possible to 
determine Mq but we can nevertheless put a lower limit 
on the ejecta mass. Models with Mo below 1O _2 5 M0 
are ruled out at the 98.3% level by the observations of 
SN 2009bb. Hence, radio observations of a relativistic 
supernovae can constrain all the physically relevant pa- 
rameters of our model. 
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7. COMPARISON WITH KNOWN SOLUTIONS 

IBlandford fc McKed (j!976f ) (hereafter BM) analyze 
the dynamics of both non-relativistic (NR) and ultra- 
relativistic (ER) blast waves in the adiabatic impulsive 
(AI) approximation as well as steady injection (SI) from 
a central power supply. A classification of the shock dy- 
namics and the corresponding synchrotron and inverse 
Compton radiation from the strong relativistic spheri- 
cal shock in an ionized ma gnetized medium is given by 
IBlandford fc McKeel (|1977l) . The dynamics of both NR 
and ER shocks are governed by the energy E of the adia- 
batic impulsive (AI) explosion; while in the case of steady 
injection of energy the dynamics is governed by Lt, where 
L is the luminosity of the central power supply. Other 
assumptions made by BM were that the shock is a strong 
one (see Footnote [2|) and that the magnetic field is not 
dynamically important. 

The external medium into which the shock wave ex- 
pands can have radial variations of density. It can ei- 
ther be of uniform density or could be stratified into 
p oc r~ k , with k = 2 for a constant velocity wind. The 
shock can be either adiabatic or radiative. In adiabatic 
shocks the radiative mechanisms are slow compared to 
hydrodynamics timescale, while in radiative shocks the 
radiative mechanisms are faster than the hydrodynamic 
timescale (measured in the observer- frame). A fully ra- 
diative blast wave may radiate away all the thermal en- 
ergy generated by the shock, if for example, the external 
medium is c o mpos ed mainly of electrons and positrons. 
ICohen e~a l. (1998) categorize shocks as "scmi-radiativc" 
in which the cooling is fast, but only a fraction of the 
energy is radiated away. This could take place in a col- 
lisionlcss shock acceleration, where the shock distributes 
the internal energy between the electrons and protons. 
After the material passes behind the shock, there is 
no coupling, effectively at the low densities behind the 
shock, between the electrons and protons. During the 
initial stage of the afterglow, the electrons may undergo 
synchrotron cooling or Compton scatter off low energy 
photons on times cales which are shorter than the dy- 
namical timesca le (lWaxmanlll997btlMeszaros et al.lll998t 
I Sari et al. 1998]), but the protons may not radiate and 
may remain hot. A narrow cooling layer may form be- 
hind the shock as in a fast cooling scenario even when the 
protons remain adiabatic and only the electrons cool. As 
the cooling parameter e (the fraction of the energy flux 
lost in the radiative layer) increases, the matter concen- 
trates in a small shell near the shock in the Newtonian 
solution; however in the ER case, even in the adiabatic 
case, the matter is concentrated in a narrow shell of width 
R/T 2 (where T = V2j, is the ultra-relativistic shock 
wave Lorentz factor in the observer's frame), but this 
concentration in a narrower denser shell increases with 
e. 

Classical fireball models of GRBs have a fully radia- 
tive stage during the (initial) 7-ray event with a radia- 
tive efficiency near unity and the implied energy is typ- 
ically E ~ 10 51-52 erg and have a bulk Lorentz factor 
of r ~ I0 2 — 1 3 . When the newly shocked electrons are 
radiating near the peak of the initial post-shock energy 
7e ~ £ e (m p /m e )r(i) (here £ e = 0(1)), they retain high 
radiat ive efficiency even afte r the GRB outburst for some 
time eszaros et al.l Q7)!)!v . If the protons establish 



and remain in cquipartition with the electrons through- 
out the entire remnant volume, the shock Lorentz factor 
evolves for a homogeneous external medium as: T cx r -3 
(jBlandford fc McKeel fl976h . On the other hand, for the 
adiabatic regime where the radiative losses do not tap 
the dominant energy reservoir in protons and magnetic 
energy but only the electrons are responsible for the ra- 
diation, one has: V oc r~ 3 l 2 for a homogeneous medium 
(|Meszaros et al.lll998l ). 

For a NR blast wave of energy E into a medium of 
density po in the Sedov- Taylor phase, it is possible to 
construct a characteristic velocity at a time t from the 
combination (E/ pot 3 ) 1 / 5 . For the relativistic problem, as 
BM point out, an additional velocity the speed of light is 
introduced into the problem. The mean energy per par- 
ticle in the shocked fluid varies as T 2 , where one factor 
of r arise from the increase in energy measured in the 
co-moving frame and the second arise from a Lorentz 
transformation into the fixed (i.e. the observer) frame. 
When most of the energy is resides in recently shocked 
particles, the total energy is proportional to T 2 R 3 . Here 
R is the current shock radius. If the total energy con- 
tained in the shocked fluid is to remain constant, then 
r 2 oc i~ 3 . One can consider a more general case in which 
the energy is supplied continuously at a rate proportional 
to a power of the time and set: 

T 2 oc t~ m for to > —I. (35) 

For the m=3 case, corresponding to an impulsive injec- 
tion of energy into the blast wave on a time scale short 
compared with R, the total energy can be shown to be 
given by: 

E = 87ru)ii 3 r 2 /17, (36) 

where w\ is the enthalpy ahead of the shock, so that with 
r 2 oc t -3 , the total energy is indeed a constant in time 
in the fixed frame. 

BM also provide a solution for an ER adiabatic blast 
wave in an external density gradient p\ oc n\ oc r~ k but 
with a cold, pressure-less external medium. This is the 
case for a blast wave propagating through a spherically 
symmetric wind. For this example they find that for k=2, 
m=l corresponds to an impulsive energy injection due to 
a blast wave in a constant velocity wind. In this case, the 
analytic solution to the adiabatic impulsive blast waves 
can be obtained for m = 3 — k > — 1, and the total energy 
in this case is: 

E = 8TT Pl T 2 t 3 /{17 - Ak). (37) 

Again this energy is constant to the lowest order in T~ 2 . 

We compare our solution to other well known analytic 
and numerical schemes in the two final Figures. The 
graphs are not really for SN 2009bb-likc systems. As 
they are for a baryon loaded explosion with a 70 of 300 
going into a wind-like CSM. There may not be any such 
spherical explosions with these kind of energies. The 
comparison in these figures are not for any physically rel- 
evant objects (or scenarios) but rather for the purpose of 
showing the consistency of our analytical solutions with 
previously reported numerical techniques and known real 
asymptotic or intermediate asymptotic solutions. 

In Figure [5] we compare our solution, in the ultra- 
relativistic regime with nearly free expansion or the 
coasting solution at early times and a BM like solution 
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at intermediate times. The blast wave solution for the 
special case of adiabatic impulsive solution for a con- 
stant velocity wind obtained in Section [3] is valid for ar- 
bitrary values of Mq and 70 and the mass loss parameter 
A. The corresponding solutions in the Blandford Mc- 
Kee analysis can be obtained from our exact solution in 
the limit of small initial mass of the high velocity ejecta 
Mq <C m(i?)7o where m{R) is the swept-up mass in the 
wind. Thus Equation (Bl) shows that our intermedi- 
ate asymptotic expressions for R and 7 have the same 
time dependence as those of the BM solution. Simi- 
larly, the opposite limit Mq 3> m(R)^/Q gives a nearly 
free expansion of the blast wave. Thus the relativistic 
CEDEX which had relatively low mass (although the 
prototype SN 2009bb was still significantly baryon loaded 
compared to classical GRBs), underwent nearly free ex- 
pa nsion initially. T his is the relati vistic analogue of 
the lChevalieil ()1982l ): iNadezhinl (|1985f ) phase of the non- 
relativistic supernovae. Eventually, an ultra-relativistic 
CEDEX would enter the Blandford-McKee phase (the 
relativistic analogue of the Sedov- Taylor phase) when it 
has swept up enough mass from the external medium 
surrounding the CEDEX. 

In Figure [6] we compare our solution, in the mildly- 
relativistic regime with a BM like solution at intermedi- 
ate times and a Snowplough phase at late times. The 
BM-like phase becomes unphysical when 7 approaches 
1 (hence (3 becomes 0) and the Snowplough phase is 
feasible only much later when (3 is less than 1. Both 
Figures also compare our solution to a direct numerical 
solution of the equations of motion computed for this 
w ork, under under the relev ant conditions as prescribed 
bv lChiang fc Dermerl (jl999f ). These Fi gures demonstrate 
that patched solutions are much worse than the solu- 
tion proposed here or even numerical solutions. While 
there can be a patch (though worse than our solution) 
between coasting and BM, there exists no patch between 
BM and Snowplough. Whil e our analysis is f or a w ind 
like {p oc r~ 2 ) medium, see iKobavashi et al.l (|1999ft for 
similar transitions seen in numerical simulations of the 
evolution of an adiabatic relativistic blastwave in an uni- 
form medium. 

Our framework may be used in the future to study the 
transition of a re lativistic bl a stwav e into a non relativis- 
tic Sedov phase. iFrail et ail ([2000D have used the Sedov 
phase for very late time calorimetry of GRB 970508. The 
relation between the blast wave radius and the observer's 
frame time also leads to the appropriate numerical factor 
connecting the two for a p cx r~ 2 exterior (See Appe ndix 
iBl). N o te the different res ults in the literature due to lSaril 
(|1997t) ; iWaxmanl (|1997aD for an uniform media. 

8. DISCUSSIONS 

Explosion dynamics where adiabatic, ultra-relativistic 
or trans-relativistic outflow takes place in a wind-like 
CSM can be classified in a parameter space spanned by 
Mq, 70 and A, i.e. the initial ejected mass, bulk Lorentz 
factor of the ejecta and the mass loss parameter of the ex- 
ternal wind established before the explosion, apart from 
the electron acceleration and magnetic field amplifica- 
tion efficiencies parametrized by e e and e#. At low 70, 
we have the non relativistic supernovae, which usually 
have considerable ejecta mass (e.g. at least 0.1 Mq, 
often much larger). If the progenitor of the SN has 



had very heavy mass loss soon before the explosion took 
place, then an initially free expansion of the shock quickly 
sweeps up a mass equal to or greater than Mq and the 
explosion enters the Sedov- Taylor like phase. Alternately 
a low ejecta mass initially would also lead to the early 
onset of the S-T phase. 

The long duration GRB afterglows are powered by ex- 
plosions that have very small initial ejecta mass, typi- 
cally 10~ 6 Mq. Comparatively, the initial mass of the 
ejecta in SN 2009bb estimated in Section |6] is consider- 
ably larger at Mq > 10 ~ 2, 5 Mq, but not as large as non- 
relativistic supernovae. These explosions are therefore 
baryon loaded. As baryon loading is a factor that affects 
the conversion of the impulsive release of energy into ki- 
netic energy of the matter around the central source, this 
can quench the emergence of the gamma-rays in a burst. 
The substantial baryon loading in SN 2009bb compared 
to classical GRBs may in fact have played a significant 
role in the circumstance that no gamma-rays were in fact 
seen from this CEDEX, despite a thorough search by a 
suite of satellites in the relevant time window. At the 
same time, the CEDEX, such as SN 2009bb, are also 
unique in that they span a region of 70 which is trans- 
relativistic. So far the LGRBs and the corresponding 
outflows were being described in the extreme relativis- 
tic limit as in BM. Here we provide a solution of the 
relativistic hydrodynamics equations which is uniquely 
tuned to the CEDEX class of objects like SN 2009bb. 
Because of the non-negligible initial ejecta mass of such 
a CEDEX, objects like this would persist in the free ex- 
pansion phase for quite a long time into their afterglow. 
Sweeping up a mass equal to that of the original ejecta 
would take considerable time, unless the mass loss scale 
in its progenitor was very intense, i.e. it had a large A. 

In this paper we also provide an analysis of the peak 
radio flux versus the product of the peak radio frequency 
and the time to rise to the radio peak. There are loci in 
this plane that are spanned by low, intermediate and high 
velocity explosions. Radio SNe, LGRBs and CEDEXs 
are seen to occupy different niches in this plane. Sim- 
ilarly, we provide expressions for the peak fluxes and 
peak radio frequencies in terms of mass-loss factors of 
the explosions. We also invert their dependence on 70 
and M in terms of the peak frequency, peak time and 
peak fluxes to interpret the parameters of the explosion. 
We have also demonstrated that a seed magnetic field 
amplified by the shock as described above plays a crucial 
role in the generation of synchrotron radiation and have 
argued (Chakraborti et. al, in prep) that the accelera- 
tion of cosmic rays to ultra-high energies is possible in 
SN 2009bb-like objects. 



We thank Alicia Sodcrberg, Abraham Locb and 
Poonam Chandra for discussions on mildly relativis- 
tic supernovae. We thank Naveen Yadav for checking 
the manuscript and Swastik Bhattacharya for help with 
mathematics. This work has made use of the Wolfram 
Mathematica computer algebra system. This research 
was supported by the TIFR 11th Five Year Plan Project 
no. 11P-409. We thank the Institute for Theory and 
Computation, Harvard University for its hospitality. We 
thank an anonymous referee for detailed comments which 
significantly extended this work, in particular for point- 



Baryon Loaded Rclativistic Blastwaves 



9 




Fig. 1. — Evolution of the blast wave radius Ri a t, determined 
from SSA fit to observed radio spectrum, as a function of the 
time t ;, s in the observer's frame. The evolution is consistent 
with nearly free expansion. Note that the observations require 
M > 10- l5 M Q . 



Fig. 3. — The flux density F v as a function of the frequency v 
and the time in the observer's frame f D (, s , as predicted from the 
model proposed in this work. The values used for the parameters, 
are those of SN 2009bb. 




Radius (cm) 



Fig. 2. — Magnetic field as a function of blast wave radius, as 
determined from SSA fits. B lue dots represen t size and magnetic 
field of radio supcrnovac from Chevalier (1998) at peak radio lumi- 
nosity. Red crosses (with 3a error-bars) give the size and magnetic 
field of SN 2009bb at different epochs, from spectral SSA fits. Red 
line gives the best B oc R~ 1 (Equation I22| l fit. 

ing out to us the Chiang and Dermer papers, which we 
were unaware of. We thank Charles Dermer for sharing 
their numerical results in these papers. We thank Tsvi 
Piran and Charles Dermer for discussions at the Annapo- 
lis GRB2010 meeting. 




Fig. 4. — Diagnostic F vv vs v v t v plot of our fiducial CEDEX in 
the nearly free expansion phase, at a distance of 40 Mpc. Solid 
lines from right to left represent loci of constant 70 = 1.005 (Red), 
1.16 (Gr een) and 3.0 (Blue), traced by keeping the LHS of Equa- 
tion [Dl] fixed at the respective values. Dashed lines from right to 
left represent loci of constant M = 0.2 (Red), 1.9 (Green) and 20 
(Blue) in un its o f 10 — 6 IV^yr" -1 , again traced by keeping the LHS 
of Equation ID2I fixed at the respective values. The cross marks 
the s et of values derived for the prototypical SN 2009bb spectrum 
fromlSodcrbcrg ct al. (2010) at 20 days, with 3<r error-bars. Radio 
observations can be used to deduce the initial bulk Lorcntz factor 
70 of a CEDEX and the pre-explosion mass loss rate M from this 
diagram. 
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Fig. 5. — Early time radial evolution of the bulk Lorcntz factor 
f(R) in a CEDEX, derived in Equation |51 of this work (red solid 
line), as compared to an initial nearly free expansion or coasting 
phase (green dotted line) and the intermediate Blandford-McKcc 
like phase (blue dashed line). Magenta squares represent a direct 
numerical sol ution of the equations of m otion following the method 
prescribed by Chiang & Dcrmcr (1999). Note the transition from 
a coasting phase to the Blandford-McKee like phase. 
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Fig. 6. — Late time radial evolution of /3(R) = v(R)/c in 
a CEDEX, derived using Equation \5\ of this work (red solid 
line), as compared to an intermediate Blandford-McKee like phase 
(blue dashed line) and a final Snowplough phase (orange dot- 
ted line). Magenta squares represent a direct numerical solution 
of the equations of mot ion following the method prescribed by 
IChiang fc D crmcr (1999]). Note the transition from the Blandford- 
McKee like phase to a Newtonian Snowplough phase. 
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APPENDIX 

A. TEMPORAL EVOLUTION OF BLASTWAVE PARAMETERS 
The evolution of the blastwave radius has already been shown to be given by 



R 



-Mo + 2Ac-y t + y/8AcM tj$ + ( M o - 2Ac l0 ty 



2A l0 



This can be substituted into into Equations §5§ to get 



70 Mo 



-M + 2Ac 70 i + ^8AcM t-y$ + (M - 2Ac 7o t) 2 
^70 



(Al) 



(A2) 



^ ^2Ac7 t + ^8AcM t^ + ( M o ~ 2Ac 70 i) 2 ^ M + 



-M + 2Ac 70 t + ^Js,AcM t~il + (M - 2Ac 7o t) 2 ^ 



47o 2 



which gives the time evolution of the Lorcntz factor 7. The amount of kinetic energy converted into thermal energy 
E is similarly obtained by substituting R(t) into Equations (J6j> to get 



E = c z 



( M - 2Ac 70 t - ,/8AcM f 7 3 + (M - 2Ac 7o t) 2 

-M + X 

^ 2 7o 



(A3) 



^ \2Acy t + ^JsAcMot^ + (M - 2Ac 70 t) 2 ^ M + 



-M + 2Ac 70 t + ^8AcM t7o + ( M o ~ 2Ac 70 t)' 
4 7 2 



Together, these three equations describe the temporal evolution of the blastwave parameters. However, for times less 
than tdec, it is sufficient to use the first few Taylor coefficients, to derive the leading order behavior of the observable 
quantities. 

B. ULTRA- RELATIVISTIC LIMIT FOR NEGLIGIBLE EJECTA MASS 



IBlan dford fc McKed (|1976f ) describe the self similar evolution of an ultra-relativistic blastwave in the 7 ^> 1 limit. 
This is the relativistic analog of the non-relativistic S e dov- T aylor solution as the ejecta mass is considered negligible 
compared to the swept up mass. IBlandford fe McKed (|1976[) provide a generalization of their solution to a blastwave 
plowing through a p oc r~ k circumstellar profile. This is of particular interest for the k = 2 constant wind profile, 
relevant for the early evolution of some GRB afterglo ws. For this situation, their solution provides a self similar 
e volution with R oc t 1 / 2 and 7 o^ " 174 fsee lPiran|[200l Equation 77 with k = 2). 

IBarenblatt fc Zel'Dovicb] (|1972|) have pointed out that self similar solutions do not merely represent specific solutions, 
but describe the intermediate asymptotic behavior of a wider class of solutions away from the boundaries of the 
independent variables. Hence, we investigate the low ejecta-mass limit (lim Mq — > 0) in our model, keeping the initial 
energy Eq = 70A/0C 2 constant and obtain from Equations (jAll and IA2[) 



I? 



2Eo 
Ac 



1/2 



,1/2 



1/4 



(Bl) 



Hence the ultra-relativistic Blandford-McKce solution fo r a p oc r~ 2 c ircumstellar media is a special limiting case of the 
general solution obtained in this work. Our Equations (jAll and IA2j) track the evolution of the blastwave from nearly 
free expansion (Mo ^ ■ m(R )jo) into the Blandford-McKee phase (Mo <C m(R)^fo). Note that unlike the self similar 
solutions (eg. Equation lBll) th ere is no divergence in quantities like the bulk Lorentz Factor at the t = boundary of 
our solution (eg. Equati on IA2[) . As the result 70 is bounded above, by total energy considerations. 

Note that Equations (|B1|) together imply an unique relation between the blastwave radius and the time in the 
observer's frame. t„h« = i?/(47 2 c), which is different from the commonly used expression (Equation [7]) given by 
iMeszaros fc Reesl (119971) b y a factor of 2. Th e corresponding relation in case of a uniform external medium has been 
argued over bv lSaril (|1997t) ; iWaxmanl (|1997al) . 

C. STIR-FRY EXPRESSIONS: RADIO SPECTRUM IN NEARLY FREE EXPANSION PHASE 



Th e relativistic outflow in the prototypical SN 2009bb was discovered from its strong radio emission ([Sodcrberg et al.l 
[2010h . Hence, more such objects may be uncovered in radio follow up of type Ibc supernovae. Equations (|27l and [28]) 



12 Chakraborti & Ray 

completely describe the early (t < td ec ) temporal evolution of the radio afterglow of a CEDEX. However, instead of 
expressing F vp and v p in CGS units in terms of many constants, it would be of use to express them in units commonly 
used by radio observ ers . Hence , substi tuting for the fundamental constants and replacing for the appropriate choice 
of c\ , C5 and Cq from iPacholczvkl (|1970( ) , we express the peak flux density as 

e B \ 9 / 14 / e e \ 5 / 7 / / x " 9/14 



(CI) 



p \O.33J V0.33/ 1 0.5/ w ; 



D 



40 Mpc 




Here, the fiducial values of the parameters are chosen from those appropriate for the prototypical SN 2009bb. Similarly, 
the temporal evolution of the SSA peak frequency is given by, 



9.5 x 



Us \ 1 f 6 B\ V U ( Ze \ 2 / U f \ - 5/14 _ nl/7 

i days J \0.33) \Q.33J {0.5) 170 j 




20 

x^ + 1 »-" 14 llTOjll&J ) GHz. (C2) 

Given the model parameters, these equations together describe the peak radio flux density of a CEDEX and the time 
evolution of its SSA peak in the nearly free expansion phase. The flux density at any frequency v and time t, can then 
be expressed as 

\ F »P ( FTtT ) if i>< v p , or 

F„(t)^{ P ) Mt) (_ (p _ 1)/2 (C3) 

in terms of the already computed F up and v p , where an electron index of p ~ 3, is appropriate for a SN 2009bb-likc 
spectrum with a optically thin spectral index of aw — 1. These expressions also be should be useful in designing radio 
surveys aimed at detecting CEDEXs. 

D. EXTRACTING BLAST WAVE PARAMETERS FROM RADIO OBSERVATIONS 

The inverse problem is that of determining the initial bulk Lorentz factor specified by 70, progenitor mass loss 
rate given by A or M and the initial ejecta mass Mo, from the radio observations. The bulk Lorentz factor may be 
determined from the radio observations using Equation Q33[) to get the simplified expression for 70 as, 

/ \ 2/19 / f \ -2/19 / , \ -2 . s9/P \ 18/19 / n \ 36/19 

.,2 _i , „«„ ( « ' (J_) ( tobs \ ( v P y 2 / F vp \ ! D \ 



75 ^ 1 + 0.225 X — — — F ' — . (Dl) 

_ \te), \0.5j \20daysy VlOGHz/ \20mJy/ Wo.Mpc/ / ' 

The result is insensitive to the cquipartition parameter a = e e /^B and filling fraction /. This may be used to reliably 

determine the initial bulk Lorentz factor of a radio detected CEDEX in the mildly relativistic, nearly free expansion 

phase (like SN 2009bb). 

An expression for A = M/v W i n d, the circumstellar density profile, set up by the mass loss from the progenitor may 
be obtained by eliminating 70 between equations and [25]) . This gives us a complicated algebraic dependence on 
F vp and v p tp. Since v p t p is a pure number 3> 1, we can expand this expression in an asymptotic series (lErdelvilll95l 
around hxafptp — > 00. This gives us the approximate expression for the mass loss rate as 




This approximate expression indicates the dependence of the inferred mass loss rate on the observational parameters, 
and makes an error of only < 10% in case of SN 2009bb, when compared to the exact expression. Given the uncer- 
tainties in the observations, we recommend the use of this expression to get an estimate of the mass loss rate from a 
CEDE X progenitor. Note that, this expression has similar scaling relations as Equati on (23) oflChevalier k, Franssonl 
(|2006[) . Hence, the mass loss rate of SN 2009bb as determined using that equation bv lSoderberg et al.l (|2010l ) remains 
approximately correct. 

The initial ejecta rest mass M cannot be estimated from radio observations in the nearly free expansion phase. 
It can only be determined when the CEDEX ejecta slows down sufficiently due to interaction with the circumstellar 
matter (Figure [TJ. Thereafter, the initial ejecta mass can be obtained from the Equation (TTJ} using the timcscale of 
slowdown tdec and the already determined A and 70. Nearly free expansion for a particular period of time, can only 
put lower limits on the ejecta mass, as shown in this work. 



Baryon Loaded Relativistic Blastwaves 



13 



E. ANGULAR SIZE EVOLUTION AND VLBI 

Very Long Base Interferometric (VLBI) measurements of the apparent angular diameters may be compared to the 
predicted value of Ri a t, as a direct test of the CEDEX model. In the nearby universe, where the luminosity distance 
and the angular distance are not significantly different, we can substitute Equation (|D1|) into Equation (|16[) to obtain 
the predicted angular diameter 9 ~ 2Ri at /D of a CEDEX as, 

c \ 1/19 /f \ "1/19 s _i/ P \9/19 / n x -1/19 



9 ~ 82X \e e ) \0.5J V 10 GHz J v 20mJyJ v 40 M P C / ^ ^ 

This expression which is insensitive to the equipartition parameter, the filling factor and even the assumed distance to 
the source, may be used in planning VLBI observations of a radio detected CEDEX and for comparing the observed 
angular sizes with those predicted from our model. 

At t bs — 81 days post explosion, the radio spectrum of SN 2009bb, constructe d from broadband Gian t Metrewave 
Radio Telescope (GMRT) and Very Large Array (VLA) observations, as given bv lSoderberg et all (|2010[ ) is well fit by 
anSSA spectrum with F vp — 10.82 ± 0.34 mJy and v v = 1.93 ± 0.07 GHz. Substituting these values into Equation 
(|E1[) we have the angular diameter as = 318 ± 12/ias, for the fiducial values of /, es and e e adop ted in this work. 
Henc e, the angular radius of 0.16 mas predicted by our model is consistent wi th the 3<7 upper l imit (Bicte nholz et aLI 



120101 ) of 0.64 mas reported from VLBI observations at around t f, s ~ 85 days. IBietenholz et al.l (|2010f ) adopt a Sedov- 



Taylor expansion for the blastwave, according to which the source will not be resolved anytime soon. However, our 
model indicates a nearly free, mildly relativistic expansion for SN 2009bb, which may soon be resolvable at the VLBI 
scale. We note, by t & s ~ 222 days, the radio spectrum had evolved to F vp = 8.3 5 ± 0.59 mJy and v v = 0.53 ± 0.04 
GHz, predicting an angular radius of around 0.51 mas according to Equation (|E1[) . We therefore predict that a 
careful VLBI observation will now successfully resolve the radio emission from SN 2009bb and confirm the nearly free 
expansion indicated by our analytic solution. Given the lower fluxes at the usually high VLBI frequencies, this may be 
a challenging observation due to inadequate sensitivity and unsuitability of self-calibration techniques at low fluxes. 
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